The coalescence remnant of spinning binary black holes 
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We compute the gravitational radiation generated in the evolution of a family of close binary black 
hole configurations, using a combination of numerical and perturbative approximation methods. We 
evolve the binaries with spins, s, aligned or counter-aligned with the orbital angular momentum from 
near the innermost stable circular orbit (ISCO) down to the final single rotating black hole. For the 
moderately spinning holes studied here the remnant Kerr black hole formed at the end of an inspiral 
process have a rotation parameter a/M ~ 0.72 + 0.32(s/m£f), suggesting it is difficult (though not 
excluded) to end up with near maximally rotating holes from such scenarios. 

PACS numbers: 04.25.Dm, 04.25. Nx, 04.30.Db, 04.70. Bw 



I. INTRODUCTION 

The coalescence of two black holes of comparable size 
will provide the most extreme dynamical tests of General 
Relativity's predictions. Such coalescences are an antici- 
pated outcome of galactic collision and core merger. The 
first galaxy with binary active galactic nuclei has, re- 
cently, been discovered by X-ray observations of NGC 
6240 with the Chandra Observatory pj. Possible evi- 
dence of merger events has also been recently presented in 
radio observations of X-shaped jet morphologies, which 
may have been produced by a sudden change in the cen- 
tral black hole's spin axis caused by a supermassive black 
hole - black hole merger 0. 

We focus on modeling the final moments of a binary 
black hole merger, which generates the last few cycles 
of radiation, carrying away a significant fraction of the 
system's energy and angular momentum. In a previous 
paperQ we have studied non-spinning binary black holes. 
Our goal here is to extend this treatment to deal with 
moderately spinning holes focusing, in particular on the 
characteristics of the final remnant black hole which the 
coalescence produces. While we generically expect col- 
lisions of binary systems with significant mass ratios, at 
a given mass systems with a mass ratio near unity are 
expected to produce the strongest gravitational radita- 
tion, producing the the largest effect on the state of the 
final black hole. We will restrict to this case here. As we 
are interested in the dynamics of these systems at their 
strongest (most non-linear) moment, we are impelled to 
apply a model which includes numerical treatment of 
Einstein's equations in their fully nonlinear form. The 
results here are thus complementary to previous studies 
in both the slow-motion weak-field post-Newtonian (PN) 
approximation, appropriate when the black holes are still 
far apart HQ and extreme mass ratio treatments based 
on near-geodesic motion^Q- 

In this paper, we adopt the viewpoint that the individ- 
ual black hole spins vary only slowly as the system ap- 
proaches the ISCO. In general then, we expect pairings of 



black holes near ISCO with arbitrary spins. We treat the 
special case of spins either aligned or anti-aligned with 
the orbital angular momentum paying particular atten- 
tion to how this addition of spin to the problem affects re- 
sulting angular momentum of the finally produced black 
hole. These configurations exclude any precessional ef- 
fects of the strong field spin interactions, but include the 
case expected to produce the most rapidly rotating rem- 
nant. We apply the Lazarus 0, H, 0, E|, approach to 
treat several systems with moderate spin, allowing us to 
calculate the radiative loss of energy and angular mo- 
mentum. 



II. THE MODEL 

The goal of the Lazarus approach is to exploit a broad 
range of analytic and numerical techniques to model 
approximately black hole binaries during the different 
stages of the coalescence. In particular, while the Close 
Limit (CL) approximation Ref. 0, 0] can describe the 
ring-down of the finally formed black hole, Numerical 
Relativity (NR) simulations are needed in order to pro- 
vide a description of the system in the strong non-linear 
merger stage of the collision. We also need a description 
for the system applicable in the Far Limit (FL) which can 
provide initial values for the NR simulations. In this pa- 
per we will apply the same CL and NR treatments which 
we used to study the nonspinning case 0, H, El EH , but 
we have extended the FL model to allow spinning black 
holes. 



A. Initial data 

We handle spinning black holes in the nearly same 
way as for the nonspinning case, with the Bowen- York- 
puncture ansatz [13|. To study plunge radiation we 
want to begin our simulations with black holes at the in- 
nermost stable quasi-circular orbit (ISCO) for spinning 
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TABLE I: Spin - ISCO data based on the effective potentiai 
method appiied to Bowen-York data 



-0.50 --0.37 --0.25 --0.12 ++0.00 



-0.08 ++0.17 



TjM 

±Y/M 

±P/M 

s/M 2 

J/M 2 

L/M 2 

MQ 

m/M 

E b /M 



7.16 

1.930 

0.230 

-0.129 

0.629 

0.887 

0.098 

0.424 

-.0159 



1.831 

0.236 

-0.0958 

0.672 

0.863 

0.105 

0.448 

-.0175 



6.29 

1.651 

0.254 

-0.0649 

0.709 

0.838 

0.118 

0.46 

-.0189 



5.70 5.05 

1.433 1.193 

0.282 0.326 

-0.0313 0.00 

0.747 0.779 



0.809 
0.136 
0.464 
-.0208 



0.779 
0.162 
0.455 
-.023 



4.70 

1.062 

0.356 

0.021 

0.799 

0.757 

0.182 

0.45 

-.0250 



4.04 

0.810 

0.451 

0.0449 

0.820 

0.730 

0.229 

0.435 

-.0279 



binary black holes as determined for the Bowen- York- 
puncture initial data. Various approaches have been de- 
veloped to compute the location and frequency of the 
ISCO. We will use results based on the effective potential 
method of 0] as generalized by Pfeiffer, Teukolsky, and 
Cook |l5j | for spinning holes. A sequence of quasi-circular 
orbit configurations determined by minimizing the bind- 
ing energy with respect to separation along sequences 
of constant mass ratio, orbital angular momentum and 
spin. The ISCO is the limit point of this sequence. We 
also note that the use of Refs. [liL IT5I ] results, strictly 
valid for the "image method", holds for the "punctures" 
since the difference outside the horizons, are very small 
in practice for the fairly detached black holes studied 
here. 
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FIG. 1: Horizon formation for S7++0.17. The solid curves 
show the numerically determined location of trapped surfaces 
including the common apparent horizon of the final black hole. 
The dashed curve indicates the independently calculated hori- 
zon location of the "background" perturbed Kerr black hole 
needed in the CL part of the treatment. 



III. 



RESULTS AND DISCUSSION 



The parameters of seven ISCO orbits are given in Ta- 
ble |U The black holes have parallel spins aligned and 
counter- alig ned with the orbital angular momentum as 
in Ref. [15|. We label the different cases by S, the z- 
component of each individual black hole's spin, scaled by 
the square of its mass. This corresponds to the scaled 
angular momentum parameter a/m for each hole. In the 
tables, I represents the proper distance between throats, 
and Et, is the binding energy for the given configuration. 
Y is the coordinate location of the punctures in confor- 
mal space (on the y-axis), P is the linear momentum 
of the holes (as measured from infinity) and are chosen 
opposite and transversal to the line joining the holes in 
conformal space, so that for our simulations the total an- 
gular momentum, J, is in the z-direction. We denote by 
s the individual spin of the holes, and L the orbital angu- 
lar momentum of the system. The angular frequency of 
the quasi-circular orbit, f2, is determined from dE^/dL. 
Finally, m denotes the individual 'puncture' (or bare) 
masses of the holes, and Ej, is the binding energy for the 
given configuration. All quantities are normalized to M, 
the total ADM mass of the system, and are computed on 
the initial time slice. 



In this section we present the results of our simulations 
beginning from the data described in the last section. As 
described in Ref. 0, we begin by numerically solving 
Einstein's equations. The difficulty of performing these 
numerical simulations limits us to brief evolutions, thus 
we next apply several criteria in determining when, or 
whether, the numerical simulation has 'linearized' and 
close-limit treatment should be applicable. We first look 
at the values of the .^-invariant |17|. which we expect to 
be within a factor of two of its background value, unity. 
We also look for independence in the waveforms on the 
time T when we transition from the numerical simula- 
tion to the close limit treatment. In particular, we ex- 
pect to find independence of the waveform phase over 
the most significant part of the waveform, as we vary T. 
We also expect to see some leveling of the dependence 
of the total radiation energy as we vary T, though this 
tends to be more difficult to achieve. For each case we 
will determine waveforms, radiation energy and angular 
momentum, and the state of the final black hole. 

The first configuration we discuss, S7++0.17, is the 
case with the strongest spins aligned with the orbital an- 
gular momentum. 

In this case we are able to run the numerical Simula- 
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TABLE II: Summary of results for spinning holes 



s 


++0.17 


++0.08 


++0.00 


--0.12 


--0.25 


T/M 


7-9 


9-11 


9-11 


10-12 


11-12 


-E/M{%) 1.7-1.9 
-j/M 2 0.06-0.08 


2.3-2.5 


2.4-2.6 


1.9-2.1 


1.9-2.1 


0.09-0.10 


0.09-0.10 


0.09-0.10 


0.09-0.10 


Mf/M 


0.982 


0.976 


0.975 


0.980 


0.980 


a f /M f 
Ji/Mf 


0.778 


0.739 


0.720 


0.679 


0.639 


0.820 


0.799 


0.779 


0.747 


0.709 


MiOpeak 


0.50-0.55 


0.45-0.50 


0.45-0.50 


0.40-0.50 


0.40-0.50 


Mujqn 


0.573 


0.552 


0.542 


0.524 


0.508 


M/tq N 


0.07714 


0.07925 


0.08006 


0.08160 


0.08285 


Tcah/M 


12 


> 15 


> 17 







tion long enough to numerically locate a common appar- 
ent horizon around the black holes. Fig. 2] shows four 
snapshots of the apparent horizon surfaces 0] on the 
orbital plane for the S7-I--I-0.17 case. The dashed line 
shows the location of the apparent horizon of the Kerr 
black hole, which we have identified as the background 
for the Close Limit calculation. The plots show that the 
two black holes are well detached at T = 0. The 'grid 
stretching' effect (due to the vanishing shift we used dur- 
ing full numerical evolution) makes them to grow in the 
coordinate space. These plots are particularly useful to 
extract qualitative information about the system. Soon 
after a common apparent horizon covers the system it 
tends to have an increasingly spherical shape. The time 
of formation of a common apparent horizon roughly gives 
an upper limit to the time at which linear theory can 
take over as has been discussed in Ref. 9] . It is expected 
that a common event horizon should have appeared some 
M 's of time earlier in the evolution of the system. On 
the independent basis of S-invariant and waveform ro- 
bustness criteria, we estimate the linearization time to 
come at T = 7 — 9M. This is consistent with the early 
formation of a apparent horizon. The key physical fea- 
ture which actually makes the close limit approximation 
effective is that the black holes share a common poten- 
tial barrier which appears even earlier than the common 
event horizon. Near the time of linearization we estimate 
the location "background horizon" at rn ~ 1.55 in the 
numerical coordinates. 

Continuing the evolution in the CL treatment, and af- 
ter two iterations of the background mass and the radi- 
ated angular momentum (starting from the initial data 
values) , we estimate the total energy lost to radiation to 
be 1.7% - 1.9% of the system's initial energy. Likewise, 
0.06 — 0.08M 2 of the initial angular momentum is lost, re- 
sulting in a final Kerr black hole formed after the plunge 
with Mf w 0.982 A/ and af/Mf w 0.778. 

The attractive nature of the spin obrit coupling makes 
the SI+- 1-0.17, with the strongest aligned component 
of spin angular momentum the highest frequency and 
closest configuration at ISCO. As we decrease s, we ex- 
pect the larger separations to require larger linearization 
times T. This expectation is consistent with our calcula- 



tions which indicate a linearization time of T = 9 — 11M 
for the 57++0.08 and SI+- 1-0.00 cases, and respectively 
T = 10 - 12M and T = 11 - 12M for the antialigned 
spin cases 57-0.12 and S7-0.25. The summarized results 
are shown in Table ITT1 Overall, energy lost to radiation 
is near 2%, but peaked at spinless case. Although the 
interaction time increases as we decrease s, suggesting 
more energy could be lost here, the mean frequency of 
the raditation decreases so that the energy is lost at a re- 
duced rate for the antialigned spin cases. The decrease in 
frequency does not affect the rate at which angular mo- 
mentum can be radiated as dramatically, which seems to 
results in a nearly flat s-dependence for the total loos 
of angular momentum. The ranges given in the table, 
are estimated from self-consistency tests based on com- 
parisons of perturbative and numerical treatments. Our 
results may also be subject to systematic errors result- 
ing from effects which may have led us to underestimate 
the linearization time. We would generally expect any 
such effect to lead to an increase amount of energy and 
angular momentum radiated, and a reduction in the fi- 
nal a/m. The results are also subject to errors in the 
initial model representing black holes at ISCO. A more 
detailed study of the non-spinning case in Q indicates 
some insensitivity in the radiation result to the precise 
location of ISCO, but further dynamical testing of this 
class of initial data should be carried out as advances in 
evolution techniques make it possible. 

A particularly interesting quantity derived from each 
of these results is the angular momentum parameter for 
the final remnant Kerr holes formed as a product of coa- 
lescence. The remnant black holes have a larger rotation 
parameter a/M, for the aligned spin cases than for the 
anti-aligned ones. In Fig. [S] error bars are estimated by 
taking the absolute maximum and minimum of observed 
the damped oscillations of the energy and angular mo- 
mentum radiated versus the transition time. This pro- 
vides us with a much larger error than self-consistency 
tests suggest, but might be more representative of the 
possible systematic errors of our approach. 

It is interesting to note here that the rotation pa- 
rameters of the final Kerr hole af/Mf, of table ILTI and 
Fig. are large, but still far from the maximally rotat- 
ing hole, suggesting it is hard to generate near maxi- 
mally rotating single holes after the plunge of two in- 
spiralling holes when the have moderate individual rota- 
tion parameters like the ones studied in this paper. A 
curve fitting to the values in the Fig. [21 gives J/M 2 ~ 
0.779+0. 2566(s/m^)-0.0941(s/m^) 2 for the initial data 
and a/M ~ 0.719 + 0.324(5/™^) for the final Kerr hole, 
where m# is the horizon mass of the individual holes 
and s its individual spin. If we extrapolate the trend 
linearly toward s/m 2 H ~ 1, the result suggests we would 
require initially aligned holes with spins > 0.85 in order 
to approach a maximally rotating remnant. Our studies 
thus enhance the conclusions based on complementary 
studies of binaries in the small mass ratio limit 7] which 
suggested that it is hard to form a maximally rotating 
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unity. Even in the equal mass case it would seem to 
take near maximally spinning black holes to produce a 
maximally spinning remnant. Of course, on a broader 
astrophysical scenario, rapidly rotating black holes may 
still be produced by accretion. 

We finally note that we also report the binding energy 
E}, of the initial configurations in table [I] for the nonspin- 
ning holes and in tabled for the binary spinning holes. 
If we interpret this binding energy as a measure of the 
energy radiated in the form of gravitational waves dur- 
ing the whole inspiral period until we reach the initial 
configuration we use to estimate the plunge radiation, 
we came to the notable conclusion that approximately as 
much energy is radiated in the few cycles following the 
plunge as in the whole prior lifetime of the binary. 



FIG. 2: Rotation parameter for the final remnant Kerr black 
hole form in coalescence of spining black holes. The plot 
shows dependence on the initial individual component of spin 
aligned with the orbital anglar momentum, and normalized to 
the horizon mass s/mjj. The initial obrital and total angular 
momenta J and Jt are also shown. 



black hole by binary merger unless the mass ratio is near 
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